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Spectroscopic evidence for the binary nature of AM CVn 
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ABSTRACT 

We analysed archival spectroscopic data of AM CVn taken with the William Herschel 
Telescope in 1996. In the literature two orbital periods for AM CVn are proposed. A 
clear S-wave in the Hei 4471, 4387 and 4143 A lines is revealed when the spectra are 
folded on the 1029 s period. No signature of this S-wave is seen when folded on 1051 s. 
Doppler tomography of the line profiles shows a clear signature of the hotspot. Using 
this we can constrain the value of K 2 to lie between 210 and 280 km s _1 . Our work 
confirms the binary nature of AM CVn beyond any doubt, establishes 1028.73 s as 
the true orbital period and supports the interpretation of AM CVn as a permanent 
superhump system. 
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1 INTRODUCTION 



AM CVn (HZ 29) was f ound as a faint blue object by 



Hu 



mason fc Zwicky (1947) and shows broad He I absorption 



lines (Greenstein & Matthews, 1957) and no hydrog en. Pe- 
riodi c brightness vari ations on a period of ~18 m in. (Smak 



|1967| ) 



and flickering (Warner & Robinson, 1972) suggested 
mass transfer in a very compact binary. A model in which a 
degenerate helium white dwarf transfers helium to another 
white dwarf, driven by the loss of angular mom entum due 



to gra vitati onal wave radiation was proposed b y Paczyhski 



(1967) and Faulkner, Flannery fc Warner (1972). See for de- 



tails on the models for AM CVn stars and a list of the 8 



currently known systems Nelemans et al. (2001) 



Since the discovery of the periodic brightness variations 
and the flickering, AM CVn has been extensively studied 
with high speed photometry (e.g |Skillman et al, 1999 ; 3ol 



heim et al., 1998) yielding multiple periodicities o n 1051 
1011 and 1029 s (see also |Solheim et al., 1991 ; Harvey 



et al., 1998). Spectroscopic observations (Patterson, Halpern 



& Shambrook, 1993) show a 13.38 hr periodicity in the skew- 
ness of the Hei absorption lines. 



According to Sofheim et al. (1998) the dominant 1051 s 
photometric period is the orbital period and the 1029 s a 
beat between the orbital period and t he 13.38 hr precession 
period of the accretion disc. However, [Skillman et al. (199!: ) 
argue that 1029 s is the orbital period and 1051 s is the 
beat period. In the latter case, the system would be similar 
to the permanent superhu mp systems among th e hydrogen 
rich cataclysmic variables (skillman et al., 1999). 

To discern between these two possibilities and to estab- 



lish beyond doubt the binary nature of AM CVn, a spec- 
troscopic signature at either period is needed. In this article 
we describe a spectroscopic study of AM CVn in which we 
found a clear signature on the 1029 s period, proving that 
this is the orbital period of the system. 



2 DATA REDUCTION 

We analysed archival spectroscopic data obtained on the 
4.2m William Herschel Telescope on February 26 and 27, 
1996, with the ISIS spectrograph. In the first night a dichroic 
was used with both the ISIS red and blue arm. Because 
of ripples in the blue part of the spectrum introduced by 
the dichroic, it was removed on the second night and only 
blue spectra were taken. Because of the ripples and the lack 
of flat-fields of high quality for the first night we mainly 
discuss the data taken on the second night. In the blue arm 
the spectra were obtained with the R600B grating covering 
3960-4760 A at a resolution of 2.0 A. 

There are 403 spectra taken on the second night (173 on 
the first), each with an integration time of 30 s. Wavelength 
calibration lamp exposures were taken approximately every 
50 spectra (about every 40 minutes). The spectra were re- 
duced using the standard data reduction package MIDAS. 
After bias subtraction and flatfielding the spectra were ex- 
tracted and wavelength calibrated. The wavelength calibra- 
tion was very stable (less than half a pixel change between 
the subsequent arc spectra) and the calibration of each arc 
spectrum was used for the 25 spectra taken before and the 
25 spectra taken after the arc without interpolation. The 
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Figure 1. Spectra of the second night averaged into 10 groups 
of 40 consecutive spectra added together after normalising the 
continuum. The spectra show the change in the shape of th e ab - 
sorption lines over the 5 hours of the observation (see Sect. 3.4). 
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Figure 2. Lomb-Scargle power spectrum of the equivalent width 
lightcurve of the 4 strongest He I lines. The peak at 3.5 cy- 
cles/hour points to an orbital period around 1030s. 



flat-field for the first night was taken at the beginning of the 
second night and had to be corrected for a small wavelength 
shift. 

The continuum of the wavelength calibrated spectra was 
fitted with a cubic spline of order 3 to line-free areas and 
used to normalise all spectra with respect to their continuum 
level. In Fig. [j] we show the resulting spectra of the second 
night after averaging them into 10 groups over the whole 
night. 



3 DATA ANALYSIS 
3.1 A period search 

In order to search for line profile variations on the uncer- 
tain orbital period of the binary, an equivalent width (EW) 
lightcurve was constructed using the combined EW of the 
four strongest He I lines. A Lomb-Scargle power spectrum 
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Figure 3. Trailed spectrogram of the He I 4387 and 4471 A line 
profiles (absorption in white). Top: after folding on the candidate 
orbital period of 1028.73 s. Bottom: after folding on 1051.2 s. 
The gray scale is chosen in order to highlight the weak emission 
components in the core of the lines. A clear emission component 
is visible in both lines when the data are folded on 1029 s only, 
indicating that this is indeed the true orbital period of AM CVn. 



(Fig. |2J) reveals a clear peak around 3.5 cycles/hour, close 
to the expected orbital period of 17 minutes. 

Because the accuracy of the candidate orbital periods 
as derived from the photometric studies of AM CVn ex- 
ceeds the accuracy of our period determination using the 
EW lightcurve, we folded the data set on the two proposed 
orbital periods, 1028.73 s (ISkillman et al., 19991) and 1051.2 s 



( polheim et al., 199c ). Both periods lie very close to the peak 
in our power spectrum. 

An emission component is seen to be moving periodi- 
cally through the line profiles when the spectra are folded 
on 1029 s, but not when they are folded on 1051 s (Fig- 
ure ^) . This S-wave is a characteristic signature in the line 
profiles of accreting binaries caused by the impact of the 
gas stream onto the accretion disc around the primar y (e.g 



Marsh, 199o| ; |Spruit fc Rutten, 1998] ; |Marsh, 1999j ). The 



phase resolved spectroscopy thus provides a definite spec- 
troscopic signature of the true orbital period of AM CVn, 
and firmly establishes its binary nature. 

The S-wave is most clearly visible in the He I lines at 
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orbital phase 



Figure 14, The radial velocities of the Hell line, derived from a 

Gaussian fit as a function of orbital phase. A sinusoidal fit is used 
to establish the absolute orbital phases and derive an upper limit 
to the radial velocity of the primary. 



4387 and 4471 A, but can also be discerned in the line at 
4143 A. The Hen lines at 4686 and 4199 A, however, do not 
show any evidence for S-waves on the orbital period. The S- 
wave is strongest moving from redshift to blueshift, but only 
barely present when moving from blue to red. This reflects 
the fact that the hotspot is mainly visible when it is at the 
front of the system as observed from Earth. Note that when 
the spot comes to the front the velocity of the stream and 
of the disk at the impact point are directed away from the 
observer and thus redshifted. 



3.2 Doppler tomography 

As a next step we applied Doppler tomography in order 
to reveal the velocity structure in the lines as well as to 
establis h the exact velocity of the S-wave. Doppler map- 



ping (Marsh & Home, 1988) uses the time dependent line 
profile shapes to reconstruct the distribution of line emis- 
sion/absorption in the corotating frame of the binary. It is 
a widely applied tool in the study of the strong emission 
lines in CVs. It has revealed a variety of structures associ- 
ated with the accretion disc around the primary in interact- 
ing binaries such as brig ht spots, spira l arms and magnetic 



streams (see for a review Marsh, 2001 ) 



So far we have used an arbitrary zero point for our or- 
bital phases since the absolute phase is not known from the 
photometry. In order to estimate the absolute orbital phase 
of the binary, where phase 0.0 is defined as superior con- 
junction of the primary star, we analysed the line profile 
behaviour when folded on the 1029 s orbital period. The 
line profiles are complex with a mix of absorption and emis- 
sion components present at any given time. The He n 4686 A 
line is an exception and consist of a single emission compo- 
nent superimposed on an absorption trough, and shows no 
evidence for an S-wave. We fitted the Hell emission compo- 
nent with a simple Gaussian, in order to measure its radial 
velocity as a function of binary phase. The radial velocity 
of the He n line shows a systematic variation with an ampli- 
tude of 53 ± 6km s _1 when fitted with a simple sine-function 
(Fig. ku). If the emission is associated with the accretion flow 



around the primary, this gives us an indication of the motion 
of the primary white dwarf. Care must be taken since asym- 
metries in the distribution of the Hell emission around the 
primary also introduces apparent radial velocity shifts. So 
rather than assuming that the derived velocity is indeed the 
projected velocity of the white dwarf (Ki), we merely use 
the relative phasing in order to construct a more reliable 
zero point for our phases. The derived velocity amplitude 
can be taken as an upper limit to the projected velocity of 
the white dwarf. Using the convention that orbital phase 
zero corresponds to superior conjunction of the primary, we 
then derive the following orbital ephemeris for AM CVn; 

T (HJD) = 2450140.6135(2) + 0.011906623(3)£ 

with the formal uncertainty of the zero point indicated be- 



tween brackets and the period taken from skillman et al 



(1999|). In all plots, the orbital phases shown are the result 



of folding using the above ephemeris. 

There are two methods of calculating Doppler tomo- 
grams from phase resolved spectroscopy. One is a straight- 
f orward back projection in conjunction with a Fourier filter 
( Home, 1991 ), the other applie s maximum entropy regular- 
isation ( |Marsh fc Home, 1988| ). Our goal is to isolate the 
weak emission components in the cores of the He I absorp- 
tion lines. To this end the filtered back projection method 
is preferred, since a maximum entropy reconstruction re- 
quires pure emission line profiles. Although more prone to 
reconstruction artifacts, the back-projection method can be 
applied to the complex line profiles of AM CVn without ad- 
ditional data processing. The results of back projecting the 
line profiles of the He I 4471 A line are plotted in Fig. ^. 
Again we compare the back projection of the line profiles 
folded on the two periods; 1029 s and 1051 s. In the case 
of 1029 s, the S-wave maps into a prominent spot in the 
Doppler tomogram that is absent when folded on 1051 s. 
The location of the emission spot is where hotspot emission 
is expected, giving us confidence in the orbital ephemeris 
that is derived and our previous interpretation of the S-wave 
as due to the hotspot on the outer edge of the accretion disc. 



3.3 The radial velocity of the mass donor 

Apart from the hotspot itself, a weak ring of emission is vis- 
ible in Fig. ^| from the rest of the accretion disc. At higher 
velocities this weak emission turns into absorption which 
produces the broad deep absorption wings in the lines. Hav- 
ing established the true orbital period of AM CVn, we can 
also constrain other system parameters using the position 
of the bright spot. The hotspot is the result of local heating 
and dissipation at the outer rim of the disc where the in- 
falling gas stream impacts. In most cases, the emission has 
the v elocity of the free falling stream (e.g. Marsh & Home 



199C). However, in other cases the velocity of the hot spot 



gas appears to be a mix of the fall velocity as well as the 



local velocity flow in the disc at the impact point (Marsh 



199C). This unfortunately means that a straightforward fit 



to the position of the hotspot using single particle trajecto- 
ries is unreliable. However, some interesting limits can still 
be obtained. First of all the velocity amplitude of the spot 
itself (403 ± 15 km s _1 ) provides an upper limit to the pos- 
sible projected velocity of the mass donor star (K2), since 
the velocity of its centre of mass cannot be larger than that 
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Figure 5. Filtered back-projection of the He I 4471 A line profiles 
on the orbital period of 1029 s as well as the photometric period 
of 1051 s. The hotspot produces a prominent spot superposed on 
weak disc emission if folded on the correct period only. (Darker 
indicates more emission). The bottom two panels are the same 
images displayed as a contour plot. Any emission significantly 
brighter than the emission from the disc ring is overplotted as a 
gray scale, revealing the strong spot in the left tomogram. 
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Figure 6. Ballistic trajectories that fit to the location of the 
hotspot. The observed position of the hotspot is plotted as a gray 
scale. Overplotted is the best fit ballistic stream trajectory, if the 
hotspot corresponds to pure free fall velocities. The corresponding 
velocities of the disc are indicated by the second trajectory above 
the ballistic one. The observed hotspot velocities need to lie in 
between these two trajectories for a given choice of K2, leading 
to possible K2 values between 210 and 265 km s — x . The Roche 
lobed shape indicates the location of the mass donor star in the 
tomogram and the three crosses the centres of mass of the two 
components and the binary system. 



of the hotspot itself. Very low values of K2 also lead to gas 
stream trajectories that are incompatible with the data. For 
a given choice of K\ , K2 , we can calculate the trajectory of 
the gas stream as well as the velocity vector of the disc flow 
in order to isolate those parts in the tomogram where hot 
spot emission is possible. We shall see later that the mass 
ratio q(= M2/M1) of AM CVn is most likely 0.087, which 
means that for a given choice of K2, K\ is derived from 
Kl = qK 2 . 

If we assume that the hotspot velocities in AM CVn 
reflect the velocity of the ballistic gas stream, we can achieve 
a good fit to the observed tomogram using K2 = 260 km s _1 
and K\ = 23 km s _1 (Fig. ^). If we relax our assumptions 
and only require that the observed velocities of the hotspot 
lie anywhere between the trajectory of the gas stream, and 
that of the disc, we can constrain K2 to lie between 210 and 
265 km s _1 . The assumed value of the mass ratio has only 
a small effect on the derived values for K2- 

The least constrained range for K2 is obtained if we 
also consider the zero point of our orbital phases to be a 
free parameter. Then, in principle, K2 can lie between 200 
and 400 km s _1 . Note, however, that in order to fit the data 
with K2 > 280 km s _1 , we need to apply an arbitrary phase 
shift in order for the stream to cross the hotspot, which also 
leads to an unrealistic impact radius that is very close to the 
white dwarf. We can thus firmly constrain the value of K2 
to lie between 210 and 280 km s _1 . 
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2450141 



Figure 7. The change in relative absorption on the redshiftcd 
versus the blueshifted side across the 13.38 hour precession cycle. 
Plotted is the line asymmetry derived from a double Gaussian fit 
as a function of HJD. A least squares sinusoidal fit is overplotted. 



3.4 Disc precession and superhumps 



Patterson et al. (1993) have found a 13.38 hr period in the 
skewness of the He I absorption lines. From Fig. |l| we see in- 
deed that the shape of the absorption lines changes during 
the 5 hrs of the total observation. This can be interpreted 
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as the varying double peaked absorption pr ofiles of a pre- 
cessing accretion disc (Patterson et al., 1993). We measured 



the skewness of the line profiles in our data by fitting multi- 
ple Gaussians to the absorption lines as follows. Since there 
is evidence for emission in the line cores, which modulates 
on the much shorter orbital period compared to the slow 
changes in the absorption lines on the precessing periods, 
the central emission was masked. Double Gaussians were 
fitted to the absorption components in order to measure the 
relative strength of the blue shifted absorption relative to 
the redshifted absorption. The spectra of both nights were 
grouped into bins of 10 consecutive spectra to improve signal 
to noise and the line asymmetry was measured accordingly. 
This double Gaussian method was found to be more reli- 
able than calculating the formal skewness of the line as a 
whole since the emission components can be masked so as 
not to distort the measurements. Fig. (?) plots the derived 
line asymmetry together with a sinusoidal least squares fit. 
We do not cover the full precession cycle, so cannot improve 
on the precession period derived by Patterson et al. (13.38 
hours) and fixed the period in the fit. The asymmetry is ex- 
pressed as the ratio between the depths of both Gaussians. 
The relative contribution of the red versus blue shifted side 
of the disc thus changes with an amplitude of 33% across the 
pr ecessing cycle. Adoptin g the very simple picture sketched 
in Patterson et al. (1993) using their (j>o as phase, we coved 
phases 0.3 - 0.75 on the first night and phases 0.2 - 0.6 on 
the second night. 



4 DISCUSSION 

The fact that AM CVn behave s very similarly to the h ydro- 
gen rich superhump systems (Patterson et al., 1993) sug- 
gests that the superhump phenomena has the same origin 
(Warner, 1995). It could be explained by precession of an 
eccentric disk, which is caused by a tidal resonance between 
the disc and the companion which is found for m ass transfer- 
ring sy stems with a mass ratio sm a ller than 0.22 Whitehurst 



(1988) up to 0.33 (Murray, 2000). Simpson fc Wood (199£ 



find a more complex behaviour of the disk, including disk 
shape changes and spiral shocks in the disk, but a similar 
'precession' period. The numerical simulations of Hirose & 



Osaki (1990, 1993) show that there is a relation between 
t he precession period, the orbital period and the mass ratio 
( |Warner, 1995] ) 



.4 



1 + q 



(1) 



where A ~ 3.73 for systems with q < 0.1. Because of the 
tidal origin of the phenomenon there is no reason to expect 
a different relation for helium discs. Using 1028.73 s for the 
orbital period and 13.38 hr for the precession period, this 
relation gives a mass ratio for AM CVn of 0.087. 

From the well known fact that for a Roche lobe filling 
star the period only depends on the mass and radius of the 
donor we can calculate the mass of the donor in AM CVn 
from the period and the mass - radius relation. From the 
two evolutionary scenarios that are proposed to lead to AM 
CVn st ars, two mass - radius relations for th e donor stars are 
found ( Warner, 1995 ; Nelemans et al., 200l|) . From the mass 
ratio the companion mass is found. In Table ^ we give the 



Table 1. System parameters for the two different mass radius 
relations: deg. for fully degenerate mass - radius relation for white 
dwarfs and semi for the semi-degener ate mass radius relation for 
the final products of heliu m CVs (from Tutukov fc Fcdorova, 1989[ 
see Nelemans et al. 2001 ). The inclination is determined assuming 



K 2 = 260 km s 



two solutio ns for AM CVn using the mass - radius relation 



M 2 


Mi 


Kl/sin i 


K2/sini 


i 


(M ) 


CM©) 


(km s- 1 ) 


(km s — ) 


o 


deg. 0.033 


0.38 


55.5 


637.4 


24.1 


semi 0.114 


1.31 


83.7 


961.8 


15.7 



Nelemans et al. (2001) 



as given m 

From the system parameters we can calculate the radial 
velocities of the two components, which we can compare to 
the values of Ki as derived from the Doppler tomograms. 
In Table |l| we give the expected values of the radial veloci- 
ties {K\,il sin i) and the resulting values for the inclination 
assuming K2 = 260 km s - . 

The derived inclinations are low, consistent with the 
fact that no eclipses are observed (which means i < 78° for 
q — 0.087). The only clue to discern the two mass - radius 
relations (and thus the formation scenario) may be the fact 
that the hotspot is obscured when moving from blue to red, 
implying not too low inclination. This would slightly favour 
the fully degenerate mass radius relation. However the mass 
ratio derived from the superhump period is too uncertain to 
draw firm conclusions. 

The value of K\ derived from the mass ratio of 0.087 
(23 km s" 1 ) seems incompatible with the interpretation of 
the radial velocity of the Hell line as being caused by the 
motion of the primary. However it could be that the mass 
ratio determined from the superhump period is incorrect for 
AM CVn. For example for a mass ratio of 0.2, the value of 
K\ would be 52 km s . This also would bring the mass of 
the primary down, which rules out the fully degenerate mass 
- radius relation, but gives somewhat more realistic value for 
the primary mass in the case of the semi-degenerate mass - 
radius relation (see for more details Nelemans et al., 2001). 

Evidence for time variability of the strength of the hot 
spot was found by folding both nights separately. Using the 
the first night only, no discernible S-wave could de identified. 
However, the second night alone, resulted in a clearly visible 
S-wave in several lines. Maybe the hot spot strength varies 
with the precession period, allthough the precession phase 
shift between the two nights is small. The poor quality of 
the data of the first night makes it for the moment impos- 
sible to draw conclusions. Better sampling of the precession 
period offers opportunities to study the strength and posi- 
tion of the hot spot as a function of the precession cycle, 
allowing a better determination of the hot spot position and 
thus the mass ratio provided adequate signal to noise can 
be obtained. 



5 CONCLUSION 

We found a clear S-wave in spectra of AM CVn when folded 
on a period of 1029 s, confirming the binary nature of AM 
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CVn beyond any possible remaining doubt and establishing 
1028.73 s as the true orbital period. 

The dominant photometrical period of 1051.2 s can be 
interpreted as the beat between the orbital period and the 
13.38 hr period that has been found in the skewness of the 
absorption lines, making the system the helium equivalent 
of the permanent superhump systems. The 13.38 hr period 
then is interpreted as the precession period of the accre- 
tion disc. The shape of the absorption lines in our data set 
changes on the same period. 

By applying Doppler tomography we were able to con- 
strain the value of K2 between 210 and 265 km s _1 , which 
given the formation channels implies a low inclination (be- 
tween ~ 15° and ~ 25°). 

Our results open a new field of studying the details of 
these intriguing ultra-compact binaries by phase resolved 
spectroscopy and Doppler tomography. 
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